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We use inelastic neutron scattering to systematically investigate the Ni-doping evolution of the 
low-energy spin excitations in BaFe2- :E Ni; C As2 spanning from underdoped antiferromagnet to over- 
doped superconductor (0.03 < x < 0.18). In the undoped state, BaFe2As2 changes from param- 
agnetic tetragonal phase to orthorhombic antiferromagnetic (AF) phase below about 138 K, where 
the low-energy (<~ 80 meV) spin waves form transversely elongated ellipses in the [H, K] plane of 
the reciprocal space. Upon Ni-doping to suppress the static AF order and induce superconductivity, 
the c-axis magnetic exchange coupling is rapidly suppressed and the momentum distribution of spin 
excitations in the [H, K] plane is enlarged in both the transverse and longitudinal directions with 
respect to the in-plane AF ordering wave vector of the parent compound. As a function of increasing 
Ni-doping x, the spin excitation widths increase linearly but with a larger rate along the transverse 
direction. These results are in general agreement with calculations of dynamic susceptibility based 
on the random phase approximation (RPA) in an itinerant electron picture. For samples near opti- 
mal superconductivity at x m 0.1, a neutron spin resonance appears in the superconducting state. 
Upon further increasing the electron-doping to decrease the superconducting transition tempera- 
ture T c , the intensity of the low-energy magnetic scattering decreases and vanishes concurrently 
with vanishing superconductivity in the overdoped side of the superconducting dome. Compar- 
ing with the low-energy spin excitations centered at commensurate AF positions for underdoped 
and optimally doped materials (x < 0.1), spin excitations in the over-doped side (x = 0.15) form 
transversely incommensurate spin excitations, consistent with the RPA calculation. Therefore, the 
itinerant electron approach provides a reasonable description to the low-energy AF spin excitations 
in BaFe2- :E Ni E As2. 



I. INTRODUCTION 

Spin excitations are thought to be a candidate for me- 
diating the electron pairing for superconductivity in un- 
conventional superconductors 1 -^. For copper oxide high- 
transition temperature (high-T c ) superconductors, super- 
conductivity arises from charge carrier doping of their an- 
tiferromagnetic (AF) Mott insulating parent compounds 
and forms a superconducting (SC) dome including under- 
doped, optimally doped, and overdoped materials 3 -. Al- 
though static AF order is suppressed in optimally doped 
superconductors, spin excitations persist throughout the 
SC dome, and vanish when superconductivity ceases to 
exist in the overdoped materials^. These results pro- 
vided compelling evidence that SC electrons in copper 
oxides are intimately associated with spin excitations 5 , 
and spin excitations may mediate electron pairing for 
superconductivity^. For iron pnictides£~— , superconduc- 
tivity can also be induced from the electron or hole dop- 
ing of their AF ordered metallic parent compounds 1 ^— . 
Because the parent compounds of iron pnictide supercon- 
ductors are metallic with hole and electron Fermi surfaces 
centered at T and M points of the reciprocal space, re- 



spectively, the AF order and superconductivity may arise 
from quasiparticle excitations between the hole and elec- 
tron pockets 1 ^ - — , much different from the local moment 
Mott physics of copper oxides 3 -. In the electron itiner- 
ant picture, the AF order in the parent compounds arises 
from Fermi surface nesting of the hole and electron pock- 
ets. Since the electron-doping that induces superconduc- 
tivity also increases the size of electron pocket near M 
points and reduces the hole-pocket size near T points, 
the static AF order is gradually suppressed with increas- 
ing electron doping and superconductivity emerges from 
the signed reversed quasiparticles excitations between the 
hole and electron pockets 1 ^—. As a consequence of open- 
ing up the electronic gaps at the hole (Ah) and electron 
(A e ) Fermi pockets in the SC state, a neutron spin reso- 
nance is expected to occur at the AF nesting wave vector 
with an energy E <2A = (\A h \ + \A e \]±2rH. Indeed, in- 
elastic neutron scattering experiments on single crystals 
of electron-doped BaFe2_ a; T E As2 (T =Co, Ni, see inset in 
Fig. 1 for the crystal structure) confirm the presence of 
the resonance^ 2 - - — . In particular, recent time-of-flight in- 
elastic neutron scattering measurements reveal that the 
high-energy (E > 100 meV) spin waves in the non-SC 
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FIG. 1: (Color online) (a) The electronic phase diagram and 
crystal structure of BaFe2- 2; Ni a ;As2. The arrows indicate 
seven doping levels studied in this experiment, (b) the co- 
linear AF structure in the FeAs-plane of BaFe2-^Ni a: As2. (c) 
and (d) Schematic pictures of const ant- energy scans along 
transverse (TR) and longitudinal (LO) directions both at 
L = and L = 1. 



BaFe.2As2 persist into the SC BaFe1.9Nio.1As2, and the 
effect of electron-doping is to form the resonance and 
modify the low-energy spin excitations 3 ^. 

If we assume that the high-energy spin excitations 
are due to localized moments and electron doping only 
affects the Fermi surface nesting and low-energy spin 
excitations^ 4 -, one can make a direct comparison between 
the measured neutron scattering wave vector profiles and 
results from the random phase approximation (RPA) cal- 
culations of the three-dimensional tight-binding model 
in the local density approximation (LDA) 35 . For ex- 
ample, the transversely elongated ellipses of the reso- 
nance in the electron-doped BaFe^-zTaAs^ -31134 from 
the spin waves in BaFe2As2 3 £ due to the enhancement 
of the intraorbital, but interband, pair scattering process 
between the d xy orbitals have been interpreted as be- 
ing more effective in giving rise to the fully gapped 5 ± - 
symmetry superconductivity 3 -. The density-functional- 
theory (DFT) calculations 3 ^ had also predicted correctly 
that the momentum anisotropy of the neutron spin reso- 
nance in the optimally hole-doped materials is rotated 
by 90° from that of the electron-doped case and be- 
comes the longitudinally elongated ellipse 38 . Moreover, 
the effect of Fermi surface nesting appears to account 
for the hole-doping evolution of the spin excitations in 
Bai_ a; K a ;Fe2As2^2.. Finally, recent neutron diffraction 
work has established, within a narrow region of x in 
BaFe2- x 2xAs2 (Fig. la), elastic transversely incommen- 
surate short-range magnetic peaks that has been hailed 
as direct evidence for the spin-density-wave order due 



to mismatch of the hole-electron pocket Fermi surface 
nesting 4 ^— . 

Given that the RPA/DFT calculations have so much 
success in describing the elastic magnetic scattering and 
low-energy spin excitations in BaFe2- a; r 2 ;As2, it is sur- 
prising that there are still no quantitative comparison 
of the doping dependence of the spin excitation pro- 
files with systematic RPA/DFT calculations. In par- 
ticular, while DFT calculations predicted that spin ex- 
citations in 7.5% electron-doped BaFe2-^T! E As2 should 
be incommensurate along the transverse direction, low- 
energy spin excitations seen by neutron scattering show 
only commensurate scattering with transversely elon- 
gated ellipses ^ 31 ' 34 . Since one can systematically carry 
out RPA/DFT calculations to obtain the imaginary part 
of the dynamic susceptibility, x'^Qi^- 1 ): as a function 
of electron-doping x in BaFe2_ 2; r :z: As2, it is important 
to compare the calculation with neutron scattering ex- 
periments focusing on wave vector and energy depen- 
dence of the low-energy spin excitations. In this arti- 
cle, we present a systematic inelastic neutron scatter- 
ing and RPA/DFT study on BaFe2- a; Ni :E As2 covering 
x = 0.03, 0.065, 0.092, 0.1, 0.12, 0.15, 0.18 shown as verti- 
cal arrows in the electronic phase diagram of Fig. 1(a) 41 . 
Consistent with earlier worke r 31 ! 34 , we find that low en- 
ergy spin excitations in BaFe2- :r Ni :r As2 are anisotropic 
and form transversely elongated ellipses at the AF or- 
der wave vector. The peak widths in both the transverse 
and longitudinal directions increase linearly with x. For 
samples near optimal superconductivity, a neutron spin 
resonance appears below T c . For samples at the over- 
doped side x = 0.15, the low-energy spin excitations form 
two transversely incommensurate peaks. Upon further 
increasing electron-doping x, the low-energy spin excita- 
tions vanish concurrently with the vanishing supercon- 
ductivity. We compare these results with RPA/DFT cal- 
culations and find them to be qualitatively similar. These 
results indicate an intimate connection between spin ex- 
citations and superconductivity, thus suggesting spin ex- 
citations play an important role for superconductivity in 
iron pnictides. 



II. EXPERIMENT 

We carried out systematic neutron scattering exper- 
iments on BaFe2- K Ni :r As2 using C5 thermal neutron 
triple-axis spectrometer at Canadian Neutron Beam Cen- 
ter in Chalker River, Canada. The final neutron energy 
was set to Ef = 14.56 meV, with pyrolytic graphite (PC) 
as monochromator, analyzer, and filters. The collima- 
tions were set to [none, 0.8°, 0.85°, 2.4°]. High qual- 
ity single crystals up to centimeter sizes were grown by 
FeAs self-flux method^ 3 -, detailed procedure and sample 
characterization were published elsewhere 4 ^. We use the 
nominal composition to represent the Ni doping level 
x. We define the wave vector Q at (q x , q y , q z ) as 
(H,K,L) = (q x a/2ir,q y b/2'iT,q z c/2iT) reciprocal lattice 
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FIG. 2: (Color online) Low energy spin excitations of the non-SC BaFe2-o;Ni I As2 (x — 0.03) with Tjv = 107 K. (a) and (b) 
Energy dependence of S(Q,(jj) at Q = (0,-1,0) and Q = (1,0,1) for T = 3,80,120 K, after subtracting the background at 
Q = (0, —1.4, 0) and Q = (1.4, 0, 1.4), respectively, (c) and (d) Dynamic spin susceptibility x"(Qi w ) after considering the Bose 
factor for L = and L = 1. (e)^(h) x"(Q,<^) of constant energy scans at T = 3, 80, 120 K and E — 7.5 meV. The solid lines 
are gaussian fits to the data. 



units (r.l.u.) using the orthorhombic unit cell, where 
a « b w 5.62 A, and c = 12.77 A. In this reciprocal space 
notation, the co-linear AF structure of the FeAs layer is 
shown in Fig. 1(b), and the in-plane AF Bragg peak po- 
sition corresponding to the Fermi surface nesting occurs 
at Q = [1,0] or [0,1] r.l.u. due to twinning. Based on 
previous work^ 1 -, we plot in Fig. 1(a) the schematic AF 
order and superconductivity schematic phase diagram. 

For the experiment, we chose seven Ni-doping levels 
spanning from the non-SC to over-doped SC samples as 
marked by the vertical arrows in Fig. 1(a). These in- 
clude x = 0.03 (lightly electron-doped non-SC sample 
with T N = 107 K)i£, 0.065 (underdoped SC sample with 
T N = 72 K and T c = 8 K ), 0.092 (nearly optimal dop- 
ing SC sample with static incommensurate short-range 
order, T N = 45 K and T c = 19 K), 0.10 (optimal dop- 
ing without AF order coexisting with superconductivity, 
T c = 20 K), 0.12 (overdoped SC sample without AF or- 
der, T c = 19 K), 0.15 (overdoped SC sample without AF 
order, T c — 14 K), 0.18 (heavily overdoped SC sample 
without AF order, T c — 9 K). In order to properly carry 
out inelastic neutron scattering experiments, we prepared 
15 to 20 pieces of single crystals for each Ni-doping lev- 
els and coaligned them using E3. The total sample mass 
of each Ni-doping level is about 10-15 grams. Similar to 
previous work^S, we aligned the samples in the [H, 0, H] 
and [0, —K, 0] scattering plane, where the c-axis is about 
23.5° from the scattering plane. In this geometry, we can 
probe the wave vector transfers both along the transverse 
and longitudinal directions near the in-plane AF posi- 
tions Q = (1,0, L) or (0,1, L) at L = 0, 1 as shown in 
Figs. 1(c) and 1(d). The co-aligned samples are loaded 
inside a top loading close cylce refrigerator for easy ex- 
change of samples. 



In order to obtain a complete picture of the doping 
evolution of the low-energy spin excitations, we also in- 
clude some results from our time-of-flight (TOF) inelastic 
neutron scattering experiments on the parent compound 
BaFe 2 As2^, optimally doped x = 0.1— and overdoped 
compounds with x — 0.15. These TOF experiments were 
carried out on the MAPS and MERLIN TOF chopper 
spectrometers at the ISIS facility, Rutherford- Appleton 
Laboratory, UKMi2£. Part of the data in the x = 0.10 
compound was collected on the TAIPAN thermal neu- 
tron triple-axis spectrometer at the Bragg Institute, Aus- 
tralian Nuclear Science and Technology Organization. 
Measurements were done with a fixed Et = 14.88 meV 
by using PC monochromator, filter and analyzer. 

To directly compare the neutron scattering results, we 
have also carried out RPA calculations based on a five- 
orbital tight-binding models. The model is obtained 
by using a LDA calculation for BaFe2- a; Ni a; As2, where 
the main effect of electron doing is to shift the chemical 
potential in a rigid band model. 



III. RESULTS 

We first describe the inelastic neutron scattering re- 
sults on BaFe2- a; Ni a ;As2 for x = 0.03. In previous work 
on the x = 0.04 with = 91 K and filamentary su- 
perconductivity, the effect of electron-doping is found to 
reduce the c-axis exchange coupling in BaFc 2 As2 and 
induce quasi- two-dimensional spin excitations 45 . The 
anisotropy spin gaps at wave vectors Q = (1,0,1) and 
(1,0,0) were 2 and 4 meV, respectively 4 - 5 . For compar- 
ison, the AF Neel temperature of the system changes 
from T N = 138 K for BaFe 2 As 2 to T N = 107 K for 
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FIG. 3: (Color online) Low energy spin excitations of the underdoped BaFe2-:cNi :E As2 (x = 0.065) with Tjv = 72 K and T c = 8 
K. (a) Energy dependence of S(Q,w) at Q = (0, -1,0) and (b) Q = (1,0, 1) for T = 3, 12,80 K. (c) and (d) Corresponding 
X"(Q, u>) for L = and L = 1. (e)~(h) of constant energy scans at E = 3 meV and (i)~(l) E — S meV for T = 3 K 

and 80 K. 



BaFei.97Nio.o3As2. In addition, there are no evidence 
for bulk superconductivity until a; = 0.05. From recent 
inelastic neutron scattering work on BaFe2As2, the spin- 
wave gaps at the Brillouin zone center and boundary are 
found to be near A(l, 0, 1) = 10 and A(l, 0, 0) = 20 meV, 
respectively-^ 6 -. Assuming that the nearest neighbors, 
next nearest neighbor, and c-axis magnetic exchange cou- 
plings are J\ a (Jib), J 2 , and J c , respectivel y 36 ' 47 , one 
can fit spin waves of BaFe2As2 using the Heisenberg 
Hamiltonian^ 6 -. The anisotropic spin gaps are related to 
the exchange couplings via 4 ^ A(l, 0, 1) = 2S[(Ji a + 2J 2 + 
J c + J s ) 2 -(J la + 2J 2 + J c ) 2 } 1 / 2 , A(1,0,0) = A(0,-1,0) = 
25[(2Ji a + 4J 2 + J s ){2Jc + Js)} 1/2 , where J s is the mag- 
netic single ion anisotropy and S is the magnetic spin 
(S = 1). For BaFe2As2, the exchange couplings are es- 
timated from global fitting of the in-plane high-energy 
spin waves using the Heisenberg Hamiltonian 36 . While 
the best fitting exchange energies are SJi a = 59.2 meV, 
SJ lb = -9.2 meV, SJ 2 = 13.6 meV, SJ C = 1.8 meV and 
S'J S =0.084 meV, these results are obtained from high- 
energy in-plane spin wave dispersions and therefore can- 
not accurately determine the c-axis exchange energy. Re- 
cently, using the in-plane magnetic exchange couplings 
determined from the high-energy spin wave data^ 6 - and 
more accurate measurements of A(l, 0, 1) and A(l, 0, 0), 
the effective c-axis exchange energy and the single ion 
anisotropy of BaFe2As2 are found to be SJ C = 0.22 meV 



and SJ S = 0.14 meV, respectively 46 . These results are 
consistent with our earlier estimations 4 ^. 

Figures 2(a) and 2(b) show the background subtracted 
energy scans 5(Q,w) at the AF zone boundary Q = 
(0,-1,0) « (1,0,0) and zone center Q = (1,0,1) for 
BaFei.97Ni .o3As2, where the background scattering was 
measured at Q = (0,-1.4,0) and Q = (1.4,0,1.4), 
respectively. The low-temperature (T = 3 K) spin- 
wave gaps at the zone center and boundary decrease to 
A(1,0,1) = 5.5 meV and A(0,-1,0) = 11 meV, re- 
spectively [Figs. 2(a) and 2(b)]. Upon warming up to 
T = 80 K ks 0.75Tat, the spin-gap values decrease rapidly 
for A(1,0,1) but remains large for A (0,-1,0). They 
completely vanish at T = 120 K « 1.127V [Fig. 2(a)]. 
Figures 2(c) and 2(d) show the energy dependence of the 
imaginary part of the dynamic susceptibility, x"(Q,cj), 
estimated via x"(Q, u) = [1 — exp(— Huj/kBT)]S(Q, lj), 
where ks is the Boltzmann constant. While the x"(Q, uS) 
shows clear spin- wave gaps at 3 K, it increases linearly 
with energy in the paramagnetic state at T = 120 K. 

From the recent work on spin excitations of opti- 
mally electron-doped BaFe1.9Nio.1As2 superconductor 34 , 
we see that the high-energy spin excitations, and there- 
fore the effective in-plane magnetic exchange energies, 
are not affected by electron-doping and superconductiv- 
ity. Assuming that the in-plane magnetic exchange cou- 
plings in BaFei.97Nio.o3As2 are unchanged from those 
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FIG. 4: (Color online) Low energy spin excitations of the nearly optimal doped BaFe2- a: Ni I As2 (x = 0.092) with Tn = 45 K 
and T c = 19 K. (a) Energy dependence of S{Q, u) at Q = (0, -1, 0) and (b) Q = (1, 0, 1) for T = 3, 25, 52 K and 155 K. (c) 
and (d) Corresponding x"(Q> w ) for L = and L = 1. (e) x"(Q,w) of constant energy scans at T = 3,25,52,80,120 K and 
E = 7.5 meV. 
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of the BaFe2As2^, we estimate SJ C = 0.066 meV and 
SJ S k meV using the newly measured A (1,0, 1) and 
A(1,0, 0) values for BaFei.97Ni .o3As 2 [Figs. 2(a) and 
2(b)]. These results are consistent with the notion that 
electron-doping in BaFe2- a; Ni a ;As2 rapidly decreases the 
out-of-plane exchange couplings 4 ^. 

To explore the in-plane momentum distribution of the 
spin excitations in BaFei.g7Nio.o3As2, we carried out 
constant-energy scans along the longitudinal [H,0,H] 
and transverse [0, — K , 0] directions at E = 7.5 meV and 
various temperatures. Figures 2(e)-2(h) show x"(Q,u;) 
along different directions at T = 3, 80, and 120 K. 



Around the AF Bragg wave vector Q = (1,0,1), spin 
waves are weakly anisotropic, being broader along the 
transverse direction than that of the longitudinal direc- 
tion [Figs. 2(e) and 2(f)]. The scattering intensity in- 
creases on warming to 80 K due to closing of the spin 
anisotropy gap. As expected, the spin excitations peak 
widths above Tn are broader than the widths of spin 
waves in the AF ordered state. For L — 0, we cannot 
find any magnetic scattering at E = 7.5 meV below Tn 
due to the presence of the spin gap. Transversely elon- 
gated paramagnetic scattering emerges at temperatures 
above T N [Figs. 2(g) and 2(h)]. 

Figure 3 summarizes the constant-energy and 
constant-Q scans for the electron underdoped 
BaFei 935M0.065AS2 with Tn = 72 K and T c — 8 
K 41 i 44 . Figures 3(a) and 3(b) show the background 
subtracted magnetic scattering in terms of energy scans 
at the AF wave vectors Q = (0, -1,0) and Q = (1,0, 1) 
at various temperatures. Although resistivity and 
diamagnetic measurements show a SC phase transition 
below T c — 8 K in this compound^ 4 -, magnetic scat- 
tering 5(Q,w) at Q = (1,0,1) shows no temperature 
dependence across T c (between 3 K and 12 K) in the 
energy range 2 meV < E < 10 meV, thus indicating 
no neutron spin resonance. Upon warming up to 80 
K (= Tn + 8 K), the paramagnetic spin excitations at 
Q = (0, —1,0) and Q = (1,0, 1) become L- independent, 
in contrast to the larger magnetic scattering intensity at 
Q = (1,0, 1) in the AF ordered state. Figures 3(c) and 
3(d) show the corresponding x"(Q,w) at Q = (0,-1,0) 
and Q = (1,0,1), respectively. In contrast to x"(Q, w ) 
for BaFei.97Nio.o3As2, the spin excitations are gapless 
at all temperatures and increase approximately linearly 
with increasing energy. 

Figures 3(e), 3(f), 3(g), and 3(h) show constant-energy 
scans below T c (3 K) and above Tn (80 K) along the lon- 
gitudinal and transverse directions at E = 3 meV and 
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L = 0, 1. At T = 3 K, there are more magnetic scatter- 
ing centered around Q = (1,0, 1) than at Q = (0, —1, 0), 
suggesting the presence of intensity modulation along the 
c-axis. On warming to 80 K above Tjv, x"(Q, w ) becomes 
similar at these two wave vectors. Figures 3(i), 3(j), 3(1), 
and 3(k) plot identical scans as those of Figs. 3(e)-3(h) 
at E = 8 meV. Similar to data at E = 3 meV, we see 
that the large differences in x"(Q,w) at L = 1 and at 
3 K essentially vanish on warming up to 80 K. These re- 
sults suggest a weak L modulation of the paramagnetic 
scattering compared with spin waves below Tjy. Com- 
paring Q-scan data at E — 3 and 8 meV, we see that 
the widths of transverse scans increase with increasing 
energy, consistent with earlier results on spin waves of 
BaFe 2 As2^. 

We now examine spin excitations of BaFe2-£ C Nl c As2 
in a narrow regime where the transverse incommen- 
surate AF order was found 4 ^. For this purpose, we 
choose BaFei.908Nio.o92As2 which has T c = 19 K and 
Tjv = 40 K 41 . Figures 4(a)-4(d) show energy depen- 
dence of the magnetic scattering S(Q,uj) and x"(Q,^) 
at Q = (0, -1, 0) and Q = (1, 0, 1) below and above T c . 
While x"(Q,u;) at Q = (0, —1,0) appears to increase ap- 
proximately linearly with increasing energy at 25 K and 
52 K, the effect of superconductivity is to suppress low- 
energy spin excitations and induce a neutron spin reso- 
nance near E = 7.5 meV [Fig. 4(c)]. At Q = (1,0,1), 
x"(Q,uj) behaves similarly except that the superconduc- 
tivity induced resonance is rather broad in energy ex- 
tending from 3 meV to 9 meV, giving a resonance energy 
of E w 6 meV [Fig. 4(d)]. These results are consistent 
with the earlier work2£. 

Since elastic scattering in BaFei.go8Nio.o92As2 forms 
static short-range transverse incommensurate AF 
order—, it is interesting to see if one can also find 
incommensurate spin excitations. Figures 4(e)-4(h) 
summarize the transverse and longitudinal scans across 
the in-plane AF wave vector at L = 1 and near the 
resonance energy of E = 7.5 meV. The intensity gain of 
the neutron spin resonance below T c is seen in both the 
transverse [Figs. 4(e) and 4(g)] and longitudinal scans 
[Figs. 4(f) and 4(h)]. Although the widths of transverse 
scans are much broader than that of the longitudinal 
scans, there is no evidence for incommensurate spin 
excitations at the resonance energy. On warming to 
higher temperatures, we see gradual reduction of the 
magnetic scattering and there are still no evidence 
for incommensurate spin excitations. In addition, the 
Q-widths of spin excitations are weakly temperature 
dependent below 120 K. Therefore, one can safely 
assume that the static short-range incommensurate 
AF order below Tjy = 40 K has little impact to spin 
excitations of BaFei.go8Nio.o92As2- 

For optimally electron-doped BaFe1.9Nio.1As2 (T c — 
20 K)2£,since recent neutron TOF experiments^ have 
already mapped out the wave vector and energy de- 
pendence of spin excitations throughout the Brillouin 
zone, we will not repeat them here but instead focus- 
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FIG. 6: (Color online) x"(Q, w ) of constant energy scans 
at T = 3, 25 K and E = 7.5 meV for the slightly overdoped 
BaFe 2 - I Ni I As 2 (x = 0.12) with T c = 19 K. 



ing on temperature dependence of the energy scans at 
Q = (1,0,1). Figures 5(a) and 5(b) show the magnetic 
scattering 5(Q,w) and x"(Q) w ) at T = 3, 25, 52, and 
159 K, respectively. While one can clearly see the pres- 
ence of a spin gap and a resonance near £«7 meV at 
3 K, the normal state x"(Q,w) is gapless and increases 
linearly with increasing energy. The x"(Q, w ) a l so de- 
creases monotonically with increasing temperature. 

Figure 6 summarizes the transverse and longitudinal 
scans around the resonance energy (E — 7.5 meV) for 
a slightly electron overdoped sample, BaFe1.8sNio.12.As2 
(T c = 19 K). This sample has the same SC transition 
temperature as that of BaFei.gosNio.o92As2 but without 
static AF order. Similar to constant-energy scans in 
the underdoped samples, we again carried out transverse 
and longitudinal scans along the [1, K, 1], [H, —1,H] and 
[H,0,H], [0, — K , 0] directions, respectively, below and 
above T c . Inspection of Figure 6 immediately reveals the 
magnetic intensity gain at wave vectors Q = (1,0, 1) and 
Q = (0,-1,0) below T c . We also note that the widths 
of transverse scans are considerably broader than that 
of the longitudinal scans. Again, there is no evidence for 
incommensurate spin excitations at the resonance energy 
to within our instrumental resolution. 

Turning our attention to a more electron overdoped 
sample BaFei.8sNio.i5As2 with T c = 14 K, we note that 
in previous polarized neutron scattering experiment on 
these samples*^, spin excitations are found to be para- 
magnetic and isotropic in space below and above T c . In 
the normal state (20 K), the x"(Q,w) at Q = (1,0, 1) is 
gapless and increases linearly with energy. Upon entering 
into the SC state (3 K), the x"(Q, w ) has a small spin gap 
of ^2 meV and a neutron spin resonance at E w 7 meV 48 . 
Similarly, the x"(Q, w) at Q — (1,0,2) has a spin gap of 
~2 meV and increases with increasing energy 4 ^. Since 
previous polarized neutron scattering experiments have 
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FIG. 7: (Color online) Low energy spin excitations of the 
overdoped BaFe 2 - a: Ni :E As2 (x = 0.15) with T c = 14 K. (a) 
~ (d) Constant energy scans at T — 3, 20 K and E = 7 
meV, where the transverse scans at T = 3 K show a small 
incommensurability with 5 — 0.098. 



already measured the energy dependence of the x"(Q> to) 
at AF wave vectors, we will not repeat them here but in- 
stead focus on the in-plane wave vector anisotropy of the 
resonance. Figures 7(a) and 7(b) show the transverse 
and longitudinal scans along the [1,K, 1] and [H,0,H] 
directions below and above T c at the resonance energy of 
E = 7 meV. While one can see an enhancement of mag- 
netic scattering below T c due to the resonance, the low- 
temperature (3 K) transverse scan also show a flat top 
consistent of having two incommensurate peaks instead 
of one Gaussian. Transverse scan along the [H,—1,H] 
direction confirm this conclusion and show two clear in- 
commensurate peaks at (—5, — 1, —S) and (5, —1,(5) with 
6 = 0.098 [Fig. 7(c)]. On the other hand, longitudinal 
scans show commensurate peaks centered at the AF wave 
vector [Figs. 7(b) and 7(d)]. 

To further probe the possible incommensurate mag- 
netic excitations, we carried out additional transverse 
scans at different energies and compare these results with 
cuts from the neutron TOF measurements on the same 
sample. Figures 8(a) and 8(b) show transverse scans 
along the [1,K, 1] direction at E = 5 and 9 meV, re- 
spectively. In the low-temperature SC state, there are 
two clear incommensurate peaks as shown by the solid 
line fits using two Gaussians. From these Gaussian fits 
to the data, we can obtain energy dependence of the 
incommensurability, giving 5 = 0.073 for E — 5 meV 
and 8 = 0.121 for E = 9 meV at 3 K. In addition, 
the incommensurate peaks appear to be more robust in 
the SC state. To confirm such a conclusion, we used 
neutron TOF chopper spectrometer MERLIN at ISIS, 
Ruthcrford-Appleton Laboratory, to measure spin exci- 
tations in BaFei.85Nio.i5As2- The incident beam energy 
was set to Ei = 25 meV, and the sample was aligned such 
that the angles between and c-axis are 9 = 7°, 11°, 16°. 
The two-dimensional wave vector dependent profile of the 
spin excitations for 6 meV < E < 8 meV at L = 1 
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FIG. 8: (Color online) (a) and (b) Transverse Q scans at 
L = 1, E = 5 meV and 9 meV below T c and above T c . 
(c) and (d) are constant energy cuts from the TOF data on 
MERLIN with energy range of 6 meV< E <8 meV and 0.7 < 
L < 1.3 at T = 5 K and 20 K. The center bar show the 
incommensurability S— 0.073, 0.121 for E — 5 meV and 9 
meV, 5= 0.102 for the TOF cut around 7 meV, and 6= 0.235 
for RPA calculation, respectively. 



can be covered in this arrangement. The wave vector 
cuts at E = 7 meV obtained by Horace software using 
the combination of three sets of TOF data are shown in 
Figs. 8(c) and 8(d). In the SC state, the transverse scan 
has a flatfish top consistent with the triple-axis data in 
Fig. 7(a). On warming to 20 K (T = 6 + T c ), the scat- 
tering shows a broad Gaussian centered around the AF 
ordering wave vector [Fig. 8(d)]. Although these re- 
sults suggest that spin excitations of BaFei.8sNio.i5As2 
change from the commensurate in the normal state to 
incommensurate in the SC state, it remains unclear if 
the commensurate-to-incommensurate transition occurs 
at T c . 

Finally, we have searched for spin excitations in heavily 
overdoped SC BaFei.82Nio.isAs2, which has no AF order 
and T c = 9 K. For the experiment, we have co-aligned 
about 8 grams of single crystals. In spite of much ef- 
forts, we have been unable to find large enough magnetic 
scattering near the AF wave vector for energies around 
8 meV. Although this does not mean that there are no 
magnetic scattering at this energy, it does suggest a dra- 
matic reduction in magnetic scattering with increasing 
Ni-doping in the overdoped regime. 

Using a RPA for a three-dimensional 5-orbital tight- 
binding model for BaFe2As2, we have carried out calcu- 
lations of the RPA spin susceptibility x"(Q, oS) for the 
normal state. This model was introduced by Graser et 
al£^ and obtained from fits of the DFT band structure 
for BaFe2As2. The RPA spin susceptibility x"(Q,w) 
is obtained from the RPA multi-orbital susceptibility 
matrix X/i P Z2^3 Z4 for orbitals li,hih and I4, which is 
related to the bare (Lindhard) susceptibility matrix 
X11 12 13 u and the Coulomb interaction matrix U through 
X RPA (Q,w) = x°(Q,w)[l - Ux {Q,u)]~ 1 . The interac- 
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FIG. 9: (Color online) Comparison between the RPA calculation (a)~(d) and experimental results (e)~(h) of the in-plane 
anisotropic spin excitations in BaFe2- a; Ni r As2. 



tion matrix U in orbital space is defined in Graser et 
al£^ and contains on-site matrix-elements for the intra- 
and inter-orbital Coulomb repulsions U and [/', and for 
the Hunds-rule coupling and pair-hopping terms J and 
J'. For this calculation, we have used spin-rotationally 
invariant parameters V — U — 2J and J' — J with 
U = 0.8 eV and J = 0.2 eV. The effect of Ni substitution 
in BaFe2- a; Ni a ;As2 is assumed to be electron-doping via 
the rigid band shift. 

The top row in Fig. 9 shows results of this calculation 
obtained for an energy E = 8 meV for different electron 
dopings. As the doping increases from (a) to (d), one 
clearly sees an enhancement of the anisotropy in spin 
excitations (transverse elongation), qualitatively similar 
to the experimental results. For a doping of x = 0.15, 
we find two transverse incommensurate peaks near the 
AF wave vector. The corresponding TOF and triple-axis 
inelastic neutron scattering measurements at E — 8 meV 
are shown in Figs. 9(e)-9(h)2ii2£. The data for x = 0.065 
were mapped out on the C5 triple-axis spectrometer with 
E = 8 meV and L=latT = 3K [Fig. 9(f)]. Figure 
9(e) is obtained by doing Q-cut of energy range 7 meV 
< E < 11 meV from the x = compound measured on 
MAPS at T = 7 K2&. Figures 9(g) and 9(h) are Q-cuts 
of the energy range 6 meV < E < 9 meV for the optimal 
doped compound x = 0.1— and the overdoped compound 
x = 0.15 measured on MERLIN at T — 5 K, respectively. 
Spin excitations form transversely elongated ellipses that 
increase with increasing electron-doping. 

On initial inspection, it appears that the RPA calcu- 
lated x"(Q,u;) spectra in Figs. 9(a)-9(d) occupy much 
larger portion of the Brillouin zone than that of the ex- 
perimental results shown in Figs. 9(e)-9(h). In addition, 
the RPA calculated x"(Q,cj) in full width half maxi- 



mum (FHWM) and incommensurability 5 are about a 
factor of 2 larger than that of the measurements. How- 
ever, when the Ni-doping dependence of the FWHM from 
the fits to transverse and longitudinal scans was plotted 
in Fig. 10(a), we see a well-defined linear dependence 
of FWHM versus x for both directions at L = and 
L = 1. These results are also consistent with the FHWM 
deduced from the TOF measurements discussed in Fig. 
9. Although the RPA calculation gives the absolute val- 
ues of FWHM that are about a factor of two larger than 
that of the experiments, the calculated spin excitation 
width also has a linear Ni-doping dependence consistent 
with the experiments. To estimate the electron-doping 
dependence of the intrinsic spin excitation Q-widths at 
E = 8 meV, we assume that spin waves in the parent 
compound BaFe2As2 are instrumental resolution limited 
at E = 8 meV— . Guassian fits to the scattering pro- 
files along the longitudinal and transverse directions at 
E = 8 meV for BaFe 2 As 2 give R LO = 0.0631 ±0.0037 and 
Rtr = 0.0811 ± 0.0053 r.l.u., respectively. These values 
are consistent with the calculated instrumental resolution 
and previous results on BaFe 2 _ 3; Co x As2^. 

Figure 10(a) shows the Ni-doping dependence of the 
longitudinal and transverse spin excitation widths at 
E = 8 meV. It is clear that the excitation widths in- 
crease linearly with increasing electron-doping along both 
directions, but with a smaller slope along the longitu- 
dinal direction. This is consistent with the electron- 
doping dependence of the x"(Q, ^) from the RPA calcu- 
lation (solid lines). Figure 10(b) plots the electron dop- 
ing dependence of the spin excitation anisotropy ratio 
A, defined as A = (Wtr — Wlo)/(Wtr + Wlo) where 
Wlo an d Wtr are intrinsic widths of spin excitations 
along the longitudinal and transverse directions, respec- 
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FIG. 10: (Color online) (a) Linear doping dependence 
of FWHM for the transverse and longitudinal wave vector 
scans, (b) Doping dependence of the in-plane anisotropy of 
spin excitations, (c) Doping dependence of peak intensity at 
Q = (0, —1,0) and (1, 0, 1) normalized by phonon intensity 
at Q = (2, 0.3, 2) and T = 80 K. The lines are guides to the 
eyes. 
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FIG. 11: (Color online) (a) and (b) Temperature dependence 
of FWHM for the transverse and longitudinal wave vector 
scans and the in-plane anisotropy for BaFe2-xNi :E As2 (x = 
0.092) at E = 7.5 meV. (c) and (d) Temperature dependence 
of in-plane anisotropy for all dopings at L = and L = 1. 



crease only slightly. Figures 11(b)- 11(d) show that the 
spin excitation anisotropy ratio A is also temperature 
independent to within the errors of our measurements. 
This is consistent with earlier neutron scattering results 
on BaFei. 85 Coo.i5As2^. 



IV. DISCUSSION AND CONCLUSIONS 



tively. When the Ni-doping increases from a; = 0.03 to 
0.15, the measured spin excitation anisotropy increases 
slightly from 0.25 to 0.35. This is consistent with the 
RPA calculation, thus suggesting that the RPA calcula- 
tion captures the essential physics in the doping depen- 
dence of the low-energy spin excitations. Figure 10(c) 
shows the Ni-doping dependence of the low-temperature 
integrated magnetic scattering near E = 8 meV, where 
we have normalized the scattering intensity for differ- 
ent doping levels via phonons at Q = (2, 0.3, 2) and 80 
K, at Q = (1,0, L) with L = 0,1. We see that the 
magnetic scattering for over-doped samples decrease sys- 
tematically with increasing doping, and appear to vanish 
around x — 0.20 near the overdoped border of super- 
conductivity dome [Fig. 1(a)]. This is consistent with 
previous nuclear magnetic resonance measurements on 
BaFe2- a; Co a: As2, where the presence of the enhanced AF 
spin excitations appears to be intimately associated with 
the dome of superconductivity, and vanish for overdoped 
sample without superconductivity^. 

To determine if the widths of spin excitations are af- 
fected by temperature, we show in Fig. 11(a) the tem- 
perature dependence of the FWHM of the longitudinal 
and transverse spin excitations for BaFei.go8Nio,o92As2 
at E = 7.5 meV. On warming from 3 K to 120 K, the 
spin excitation widths are essentially unchanged, and in- 



Since the discovery of the static AF order in the parent 
compounds of iron-based superconductors^— , a central 
question has been whether the static AF order and asso- 
ciated spin excitations can be entirely described by Fermi 
surface nesting between the hole pockets near V point and 
electron pockets M points of the Brillouin zon o 50 ' 51 , or 
requires local moments as in the case of copper oxide 
superconductors^— . From recent neutron scattering 
experiments on spin waves in BaFe2As2 3 ™ and spin ex- 
citations in optimally electron doped BaFei.gNio.iAs2^, 
we see that the effect of electron-doping on BaFe2As2 
modifies spin waves below 100 meV, and does not change 
high-energy spin excitations. These results suggest that 
spin excitations in BaFe2-xNi x As2 have both itinerant 
and localized components^-. Given the general agree- 
ment on the evolution of the observed spin excitation 
anisotropy and the RPA calculations based on rigid band 
shift and itinerant electrons, one may assume that itiner- 
ant electrons and Fermi surface nesting only affect low- 
energy spin excitations, and high-energy spin waves and 
excitations arise mostly from local moments and elec- 
tron correlations. These results are consistent with the 
idea that the transversely elongated spin excitations in 
BaFe2_ x Ni x As2 are mostly due to intraorbital, but in- 
terband, scattering processes in cases without perfect 
nesting^! 

The experimental observation of a rapid reduction in 
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the spin wave anisotropy gap upon electron-doping in 
BaFe2_ a ;Ni :r AS2 confirms the earlier results that the dom- 
inate effect of electron doping is to reduce c-axis spin ex- 
change coupling and change three-dimensional spin waves 
into quasi-two-dimensional spin excitations^. Based on 
our systematic measurements of the transverse and longi- 
tudinal widths of spin excitations for different Ni-doping 
levels x, we find that the intrinsic excitation widths in- 
crease linearly with x. This is consistent with the RPA 
calculation assuming that the effect of Ni-doping is to 
increase electron Fermi pocket size via rigid band shift. 
Although the RPA calculated spin susceptibility occupies 
a larger part of the Brillouin zone than that of the ex- 
periments, their electron-doping dependences are rather 
similar. These results suggest that itinerant electron pic- 
ture and Fermi surface nesting can capture an important 
part of the physics in these materials. Similar RPA calcu- 
lations for much higher energy spin excitations (E > 100 
meV) give results that disagree with our measurements, 
thus confirming the notion that the high-energy spin ex- 
citations in BaFe2- a: Ni :c As2 may originate from the local 
moments instead of Fermi surface nesting and itinerant 
electrons^. 

In conclusion, we use inelastic neutron scattering 
to demonstrate the presence of anisotropic in-plane 
spin excitations at low energies in electron doped 
BaFe2-xNi x As2. The excitation widths in both the 
transverse and longitudinal directions increase linearly 
with doping level x, and having slightly larger slope in the 
transverse direction. In the overdoped side with x — 0.15, 



we find evidence for the low-energy transverse incommen- 
surate spin excitations consistent with the RPA calcula- 
tion. Therefore, the in-plane spin excitation anisotropy 
increases slightly with doping. For samples near opti- 
mal superconductivity, a neutron spin resonance appears 
below T c . However, the intensity of spin excitations de- 
creases with increasing doping for samples beyond opti- 
mal superconductivity, and vanishes near the overdoped 
border of the SC dome. Therefore, our data support 
the view that low-energy spin excitations are controlled 
by Fermi surface nesting and itinerant electrons, and are 
important for superconductivity of iron pnictides. 
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